Machine Learning Deciphered Molecular Mechanistics with Accurate Kinetic and Thermodynamic Prediction.

Time-lagged independent component analysis (tICA) and the Markov state model (MSM) have been extensively employed for extracting conformational dynamics and kinetic community networks from unbiased trajectory ensembles. However, these techniques may not be the optimal choice for elucidating transition mechanisms within low-dimensional representations, especially for intricate biosystems. Unraveling the association mechanism in such complex systems always necessitates permutations of several essential independent components or collective variables, a process that is inherently obscure and may require empirical knowledge for selection. To address these challenges, we have implemented an integrated unsupervised dimension reduction model: uniform manifold approximation and projection (UMAP) with hierarchy density-based spatial clustering of applications with noise (HDBSCAN). This approach effectively generates low-dimensional configurational embeddings. The hierarchical application of this architecture, in conjunction with MSM, reveals global kinetic connectivity while identifying local conformational states. Consequently, our methodology establishes a multiscale mechanistic elucidation framework. Leveraging the benefits of the uniform sample distribution and a denoising approach, our model demonstrates robustness in preserving global and local data structures compared to traditional dimension reduction methods in the field of MD analysis area. The interpretability of hyperparameter selection and compatibility with downstream tasks are cross-validated across various simulation data sets, utilizing both computational evaluation metrics and experimental kinetic observables. Furthermore, the predicted Mcl1-BH3 association kinetics (0.76 s-1) is in close agreement with surface plasmon resonance experiments (0.12 s-1), affirming the plausibility of the identified pathway composed of representative conformations. We anticipate that the devised workflow will serve as a foundational framework for studying recognition patterns in complex biological systems. Its contributions extend to the exploration of protein functional dynamics and rational drug design, offering a potent avenue for advancing research in these domains.

Computational drug development for membrane protein targets.

The application of computational biology in drug development for membrane protein targets has experienced a boost from recent developments in deep learning-driven structure prediction, increased speed and resolution of structure elucidation, machine learning structure-based design and the evaluation of big data. Recent protein structure predictions based on machine learning tools have delivered surprisingly reliable results for water-soluble and membrane proteins but have limitations for development of drugs that target membrane proteins. Structural transitions of membrane proteins have a central role during transmembrane signaling and are often influenced by therapeutic compounds. Resolving the structural and functional basis of dynamic transmembrane signaling networks, especially within the native membrane or cellular environment, remains a central challenge for drug development. Tackling this challenge will require an interplay between experimental and computational tools, such as super-resolution optical microscopy for quantification of the molecular interactions of cellular signaling networks and their modulation by potential drugs, cryo-electron microscopy for determination of the structural transitions of proteins in native cell membranes and entire cells, and computational tools for data analysis and prediction of the structure and function of cellular signaling networks, as well as generation of promising drug candidates.

Structural heterogeneity of the ion and lipid channel TMEM16F.

Transmembrane protein 16 F (TMEM16F) is a Ca2+-activated homodimer which functions as an ion channel and a phospholipid scramblase. Despite the availability of several TMEM16F cryogenic electron microscopy (cryo-EM) structures, the mechanism of activation and substrate translocation remains controversial, possibly due to restrictions in the accessible protein conformational space. In this study, we use atomic force microscopy under physiological conditions to reveal a range of structurally and mechanically diverse TMEM16F assemblies, characterized by variable inter-subunit dimerization interfaces and protomer orientations, which have escaped prior cryo-EM studies. Furthermore, we find that Ca2+-induced activation is associated to stepwise changes in the pore region that affect the mechanical properties of transmembrane helices TM3, TM4 and TM6. Our direct observation of membrane remodelling in response to Ca2+ binding along with additional electrophysiological analysis, relate this structural multiplicity of TMEM16F to lipid and ion permeation processes. These results thus demonstrate how conformational heterogeneity of TMEM16F directly contributes to its diverse physiological functions.

Cryo-EM structures of ClC-2 chloride channel reveal the blocking mechanism of its specific inhibitor AK-42.

ClC-2 transports chloride ions across plasma membranes and plays critical roles in cellular homeostasis. Its dysfunction is involved in diseases including leukodystrophy and primary aldosteronism. AK-42 was recently reported as a specific inhibitor of ClC-2. However, experimental structures are still missing to decipher its inhibition mechanism. Here, we present cryo-EM structures of apo ClC-2 and its complex with AK-42, both at 3.5 Å resolution. Residues S162, E205 and Y553 are involved in chloride binding and contribute to the ion selectivity. The side-chain of the gating glutamate E205 occupies the putative central chloride-binding site, indicating that our structure represents a closed state. Structural analysis, molecular dynamics and electrophysiological recordings identify key residues to interact with AK-42. Several AK-42 interacting residues are present in ClC-2 but not in other ClCs, providing a possible explanation for AK-42 specificity. Taken together, our results experimentally reveal the potential inhibition mechanism of ClC-2 inhibitor AK-42.

Applications and prospects of cryo-EM in drug discovery.

Drug discovery is a crucial part of human healthcare and has dramatically benefited human lifespan and life quality in recent centuries, however, it is usually time- and effort-consuming. Structural biology has been demonstrated as a powerful tool to accelerate drug development. Among different techniques, cryo-electron microscopy (cryo-EM) is emerging as the mainstream of structure determination of biomacromolecules in the past decade and has received increasing attention from the pharmaceutical industry. Although cryo-EM still has limitations in resolution, speed and throughput, a growing number of innovative drugs are being developed with the help of cryo-EM. Here, we aim to provide an overview of how cryo-EM techniques are applied to facilitate drug discovery. The development and typical workflow of cryo-EM technique will be briefly introduced, followed by its specific applications in structure-based drug design, fragment-based drug discovery, proteolysis targeting chimeras, antibody drug development and drug repurposing. Besides cryo-EM, drug discovery innovation usually involves other state-of-the-art techniques such as artificial intelligence (AI), which is increasingly active in diverse areas. The combination of cryo-EM and AI provides an opportunity to minimize limitations of cryo-EM such as automation, throughput and interpretation of medium-resolution maps, and tends to be the new direction of future development of cryo-EM. The rapid development of cryo-EM will make it as an indispensable part of modern drug discovery.

Cryo-EM structure of human heptameric pannexin 2 channel.

Pannexin 2 (Panx2) is a large-pore ATP-permeable channel with critical roles in various physiological processes, such as the inflammatory response, energy production and apoptosis. Its dysfunction is related to numerous pathological conditions including ischemic brain injury, glioma and glioblastoma multiforme. However, the working mechanism of Panx2 remains unclear. Here, we present the cryo-electron microscopy structure of human Panx2 at a resolution of 3.4 Å. Panx2 structure assembles as a heptamer, forming an exceptionally wide channel pore across the transmembrane and intracellular domains, which is compatible with ATP permeation. Comparing Panx2 with Panx1 structures in different states reveals that the Panx2 structure corresponds to an open channel state. A ring of seven arginine residues located at the extracellular entrance forms the narrowest site of the channel, which serves as the critical molecular filter controlling the permeation of substrate molecules. This is further verified by molecular dynamics simulations and ATP release assays. Our studies reveal the architecture of the Panx2 channel and provide insights into the molecular mechanism of its channel gating.

Molecular basis of ligand selectivity for melatonin receptors.

Sleep disorders in adults are related to adverse health effects such as reduced quality of life and increased mortality. About 30-40% of adults are suffering from different sleep disorders. The human melatonin receptors (MT1 and MT2) are family A G protein-coupled receptors that respond to the neurohormone melatonin MEL which regulates circadian rhythm and sleep. Many efforts have been made to develop drugs targeting melatonin receptors to treat insomnia, circadian rhythm disorders, and even cancer. However, designing subtype-selective melatonergic drugs remains challenging due to their high similarities in both sequences and structures. MEL (a function-selective compound with a bulky ß-naphthyl group) behaves as an MT2-selective antagonist, whereas it is an agonist of MT1. Here, molecular dynamics simulations were used to investigate the ligand selectivity of MT receptors at the atomic level. We found that the binding conformation of MEL differs in different melatonin receptors. In MT1, the naphthalene ring of MEL forms a structure perpendicular to the membrane surface. In contrast, there is a 130° angle between the naphthalene ring of MEL and the membrane surface in MT2. Because of this conformational difference, the MEL leads to a constant water channel in MT1 which activates the receptor. However, MEL hinders the formation of continuous water channels, resulting in an inactive state of MT2. Furthermore, we found that A1173.29 in MT2 is a crucial amino acid capable of hindering the conformational flip of the MEL molecule. These results, coupled with previous functional data, reveal that although MT1 and MT2 share highly similar orthosteric ligand-binding pockets, they also display distinctive features that could be used to design selective compounds. Our findings provide new insights into functionally selective melatonergic drug development for sleep disorders.

Asymmetric opening of the homopentameric 5-HT3A serotonin receptor in lipid bilayers.

Pentameric ligand-gated ion channels (pLGICs) of the Cys-loop receptor family are key players in fast signal transduction throughout the nervous system. They have been shown to be modulated by the lipid environment, however the underlying mechanism is not well understood. We report three structures of the Cys-loop 5-HT3A serotonin receptor (5HT3R) reconstituted into saposin-based lipid bilayer discs: a symmetric and an asymmetric apo state, and an asymmetric agonist-bound state. In comparison to previously published 5HT3R conformations in detergent, the lipid bilayer stabilises the receptor in a more tightly packed, 'coupled' state, involving a cluster of highly conserved residues. In consequence, the agonist-bound receptor conformation adopts a wide-open pore capable of conducting sodium ions in unbiased molecular dynamics (MD) simulations. Taken together, we provide a structural basis for the modulation of 5HT3R by the membrane environment, and a model for asymmetric activation of the receptor.

International Union of Basic and Clinical Pharmacology. CX. Classification of Receptors for 5-hydroxytryptamine; Pharmacology and Function.

5-HT receptors expressed throughout the human body are targets for established therapeutics and various drugs in development. Their diversity of structure and function reflects the important role 5-HT receptors play in physiologic and pathophysiological processes. The present review offers a framework for the official receptor nomenclature and a detailed understanding of each of the 14 5-HT receptor subtypes, their roles in the systems of the body, and, where appropriate, the (potential) utility of therapeutics targeting these receptors. SIGNIFICANCE STATEMENT: This review provides a comprehensive account of the classification and function of 5-hydroxytryptamine receptors, including how they are targeted for therapeutic benefit.

The G protein coupled receptor CXCR4 designed by the QTY code becomes more hydrophilic and retains cell signaling activity.

G protein-coupled receptors (GPCRs) are vital for diverse biological functions, including vision, smell, and aging. They are involved in a wide range of diseases, and are among the most important targets of medicinal drugs. Tools that facilitate GPCR studies or GPCR-based technologies or therapies are thus critical to develop. Here we report using our QTY (glutamine, threonine, tyrosine) code to systematically replace 29 membrane-facing leucine, isoleucine, valine, and phenylalanine residues in the transmembrane α-helices of the GPCR CXCR4. This variant, CXCR4QTY29, became more hydrophilic, while retaining the ability to bind its ligand CXCL12. When transfected into HEK293 cells, it inserted into the cell membrane, and initiated cellular signaling. This QTY code has the potential to improve GPCR and membrane protein studies by making it possible to design functional hydrophilic receptors. This tool can be applied to diverse α-helical membrane proteins, and may aid in the development of other applications, including clinical therapies.

Enhancing the Signaling of GPCRs via Orthosteric Ions.

G protein-coupled receptors play essential roles in cellular processes such as neuronal signaling, vision, olfaction, tasting, and metabolism. As GPCRs are the most important drug targets, understanding their interactions with ligands is of utmost importance for discovering related new medicines. In many GPCRs, an allosteric sodium ion next to the highly conserved residue D2.50 has been proposed to stabilize the inactive receptor state by mediating interactions between transmembrane helices. Here, we probed the existence of internal and functionally important sodium ions in the dopamine D2 receptor, using molecular dynamics simulations. Besides a new sodium ion at the allosteric ligand binding site, we discovered an additional sodium ion, located close to the orthosteric ligand binding site. Through cell-based activation assays, the signaling of D2 receptor with site-specific mutations was tested against a series of chemically modified agonists. We concluded an important structural role of this newly discovered orthosteric sodium ion in modulating the receptor signaling: It enables the coordination of a polar residue in the ligand binding site with an appropriately designed agonist molecule. An identical interaction was also observed in a recently released high-resolution crystal structure of mu-opioid receptor, which was reresolved in this work. Probably because of similar interactions, various metal ions have been found to increase the signaling of many other GPCRs. This unique principle and strategy could be used to optimize the drug activity of GPCR. Our findings open a new mechanistic opportunity of GPCR signaling and help design the next generation of drugs targeting GPCRs.

Advancing Drug Discovery via Artificial Intelligence.

Drug discovery and development are among the most important translational science activities that contribute to human health and wellbeing. However, the development of a new drug is a very complex, expensive, and long process which typically costs 2.6 billion USD and takes 12 years on average. How to decrease the costs and speed up new drug discovery has become a challenging and urgent question in industry. Artificial intelligence (AI) combined with new experimental technologies is expected to make the hunt for new pharmaceuticals quicker, cheaper, and more effective. We discuss here emerging applications of AI to improve the drug discovery process.

Computational modeling of the olfactory receptor Olfr73 suggests a molecular basis for low potency of olfactory receptor-activating compounds.

The mammalian olfactory system uses hundreds of specialized G-protein-coupled olfactory receptors (ORs) to discriminate a nearly unlimited number of odorants. Cognate agonists of most ORs have not yet been identified and potential non-olfactory processes mediated by ORs are unknown. Here, we used molecular modeling, fingerprint interaction analysis and molecular dynamics simulations to show that the binding pocket of the prototypical olfactory receptor Olfr73 is smaller, but more flexible, than binding pockets of typical non-olfactory G-protein-coupled receptors. We extended our modeling to virtual screening of a library of 1.6 million compounds against Olfr73. Our screen predicted 25 Olfr73 agonists beyond traditional odorants, of which 17 compounds, some with therapeutic potential, were validated in cell-based assays. Our modeling suggests a molecular basis for reduced interaction contacts between an odorant and its OR and thus the typical low potency of OR-activating compounds. These results provide a proof-of-principle for identifying novel therapeutic OR agonists.

New Binding Sites, New Opportunities for GPCR Drug Discovery.

Many central biological events rely on protein-ligand interactions. The identification and characterization of protein-binding sites for ligands are crucial for the understanding of functions of both endogenous ligands and synthetic drug molecules. G protein-coupled receptors (GPCRs) typically detect extracellular signal molecules on the cell surface and transfer these chemical signals across the membrane, inducing downstream cellular responses via G proteins or ß-arrestin. GPCRs mediate many central physiological processes, making them important targets for modern drug discovery. Here, we focus on the most recent breakthroughs in finding new binding sites and binding modes of GPCRs and their potentials for the development of new medicines.

Exploring a new ligand binding site of G protein-coupled receptors.

Identifying a target ligand binding site is an important step for structure-based rational drug design as shown here for G protein-coupled receptors (GPCRs), which are among the most popular drug targets. We applied long-time scale molecular dynamics simulations, coupled with mutagenesis studies, to two prototypical GPCRs, the M3 and M4 muscarinic acetylcholine receptors. Our results indicate that unlike synthetic antagonists, which bind to the classic orthosteric site, the endogenous agonist acetylcholine is able to diffuse into a much deeper binding pocket. We also discovered that the most recently resolved crystal structure of the LTB4 receptor comprised a bound inverse agonist, which extended its benzamidine moiety to the same binding pocket discovered in this work. Analysis on all resolved GPCR crystal structures indicated that this new pocket could exist in most receptors. Our findings provide new opportunities for GPCR drug discovery.

Single-Vesicle Assays Using Liposomes and Cell-Derived Vesicles: From Modeling Complex Membrane Processes to Synthetic Biology and Biomedical Applications.

The plasma membrane is of central importance for defining the closed volume of cells in contradistinction to the extracellular environment. The plasma membrane not only serves as a boundary, but it also mediates the exchange of physical and chemical information between the cell and its environment in order to maintain intra- and intercellular functions. Artificial lipid- and cell-derived membrane vesicles have been used as closed-volume containers, representing the simplest cell model systems to study transmembrane processes and intracellular biochemistry. Classical examples are studies of membrane translocation processes in plasma membrane vesicles and proteoliposomes mediated by transport proteins and ion channels. Liposomes and native membrane vesicles are widely used as model membranes for investigating the binding and bilayer insertion of proteins, the structure and function of membrane proteins, the intramembrane composition and distribution of lipids and proteins, and the intermembrane interactions during exo- and endocytosis. In addition, natural cell-released microvesicles have gained importance for early detection of diseases and for their use as nanoreactors and minimal protocells. Yet, in most studies, ensembles of vesicles have been employed. More recently, new micro- and nanotechnological tools as well as novel developments in both optical and electron microscopy have allowed the isolation and investigation of individual (sub)micrometer-sized vesicles. Such single-vesicle experiments have revealed large heterogeneities in the structure and function of membrane components of single vesicles, which were hidden in ensemble studies. These results have opened enormous possibilities for bioanalysis and biotechnological applications involving unprecedented miniaturization at the nanometer and attoliter range. This review will cover important developments toward single-vesicle analysis and the central discoveries made in this exciting field of research.

Expression, Biochemistry, and Stabilization with Camel Antibodies of Membrane Proteins: Case Study of the Mouse 5-HT3 Receptor.

There is growing interest in the use of mammalian protein expression systems, and in the use of antibody-derived chaperones, for structural studies. Here, we describe protocols ranging from the production of recombinant membrane proteins in stable inducible cell lines to biophysical characterization of purified membrane proteins in complex with llama antibody domains. These protocols were used to solve the structure of the mouse 5-HT3 serotonin receptor but are of broad applicability for crystallization or cryo-electron microscopy projects.

Two-dimensional crystallization of the mouse serotonin 5-HT3A receptor.

The mouse serotonin 5-HT3A receptor is a homo-pentameric ligand-gated ion channel (pLGIC) mediating fast excitatory neurotransmission in the central nervous system. The molecular mechanism of ion permeation of 5-HT3A receptors triggered by the neurotransmitter serotonin is not yet fully understood. The recent X-ray structure of the mouse serotonin 5-HT3A receptor in complex with a stabilizing nanobody revealed for the first time the entire structure of a mammalian pLGIC in detergent. Structural information of the receptor in a lipid bilayer however is still limited primarily due to the lack of 2D crystals of the receptor in a lipid bilayer. Here we present our results on the formation and improvement of diffracting 2D crystals of the mouse 5-HT3A by limited proteolysis and addition of conformational nanobodies.